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ABSTRACT 

Compact Galactic binaries where at least one member is a white dwarf or neutron star constitute 
the majority of individually detectable sources for future low-frequency space-based gravitational- wave 
(GW) observatories; they also form an unresolved continuum, the dominant Galactic foreground at 
frequencies below a few mHz. Due to the paucity of electromagnetic observations, the majority of 
studies of Galactic-binary populations so far have been based on population-synthesis simulations. 
However, recent surveys have reported several new detections of white-dwarf binaries, providing new 
constraints for population estimates. In this article, we evaluate the impact of revised local densities 
of interacting white-dwarf binaries on future GW observations. Specifically: we consider five scenarios 
that explain these densities with different assumptions on the formation of interacting systems; we 
simulate corresponding populations of detached and interacting white-dwarf binaries; we estimate the 
number of individually detectable GW sources and the magnitude of the confusion-noise foreground, as 
observed by space-based detectors with 5- and 1-Mkm arms. We confirm earlier estimates of thousands 
of detached-binary detections, but project only few ten to few hundred detections of interacting 
systems. This reduction is partly due to our assessment of detection prospects, based on the iterative 
identification and subtraction of bright sources with respect to both instrument and confusion noise. 
We also confirm earlier estimates for the confusion-noise foreground, except in one scenario that 
explains smaller local densities of interacting systems with smaller numbers of progenitor detached 
systems. 

Subject headings: gravitational waves, Galactic binaries, AM CVn, data analysis, LISA 



1. INTRODUCTION 



Compact Galactic binaries with periods shorter than 
approximately five hours are of considerable interest as 



gravi tational wave (GW) sources ( Nelemans|2009 Marsh 
20111. Such binaries include detached double white 



dwarf (DDWD) systems, short-period cataclysmic vari- 
ables (CVs), and AM CVn systems (so nam ed after 
their prototype, AM Canum Venaticorum; see |Solheim | 
20101. These last are short-period binaries in which a 
white dwarf (WD) is accreting He-rich material from an 
H-deficient companion. AM CVn systems are defined 
spectroscopically by the presence of strong absorption 
or emission He lines and the absence of H lines, which 
instead characterize "standard" CVs. In this paper we 
focus on those binaries that are evolving toward, or have 
evolved from, extremely short orbital periods, on the or- 
der of minutes. These include DDWD and AM CVn 
systems, sometimes referred to as "ultra-compact" bi- 
nary systems. Such binaries have received considerable 
attention, among other reasons, as possible proge nitors 
for Type la supernovae, R Coronae Borealis stars (Wcb- 
bhjk]|l984i JIben fe Tut ukov|J^4| , and so-called ^Ia"" 
explosions ( jBildsten et al.||2007p . Our understanding of 
their formation and evolution remains poor. 
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Together with ultra-compact X-ray binaries, DDWD 
and AM CVn systems will be the most numerous and 
individually brightest GW sources in the Galaxy for 
low-frequency (10 -0.1 Hz) space-based GW detec- 
tors. Within this frequency band, observations suggest 
0.5-2.5 x 10 7 DDWD systems and theoretical e stimates 
predict 10 5 -10 7 AM CVn systems in our Galax y ( Maxted 
fc Marsh|[l999| |Nelemans et al.||2001a|c||2004[ ). 



1'he aim ot this paper is to characterize the detectable 
GW emission from compact binaries in light of recent 
population estimates based on new observations. The 
first estimates were based solely on population-synthesis 
simulations; see App. [X] for a detailed discussion and 
review of the literature. Recent observations of AM CVn 
systems have provided new information that constrains 
earlier estimates. For instance, the identification of AM 
CVn systems among candidates selected from the Sloan 
Digital Sky Survey ( SPSS) spectroscopic catalog allowed 



Roelofs et al. ( 2007b ) to estimate the density of AM CVn 
systems by modeling the completeness of the catalog. A 
follow-up spectroscopic survey of a color-selected subset 
of the SDSS photometric database seem s to find even 
fewer AM CVn syste ms than expected ( [Roelofs et al. 



2007b Rau et al. 2010 ), leading to local density estimates 
smaller by 10-20 than suggested by population-synthesis 
studies. 

To date, 27 AM CVn systems have been presented 
in the literature; the most recent detections have been 
achieved in synoptic s urveys, including a sho rt-period 



system from Kepler (Fontaine et al. 2011), and a 



medium-period, o utbursting system fr om the Palomar 
Transient Factory (Levitan et al. 2011). 
Similarly, surveys of low-mass WD candidates selected 
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from SDSS has led to the discovery of several short- 
period DDWDs that could be progenitors of AM CVn 



systems (iBadenes et al.1200 9; Mullally et al.|2009 

' ' ]2on\ . 



rni & van Kerkwijk||2(JlU 
ample 



Marsh et al. 
the ELM ( E xtremcly "Xow Mass] survey (Brown 
et al.|2010||Kilic et al.|2011a||Brown et al.|2012| searches 
for companions around extremely low mass WDs (~ 
0.2 Mo). The ELM survey has approximately quintu- 
pled (with 19 new systems) the known sample of DDWDs 
that will merge wit hin a Hubble time ( Kilic et al.^Ollb" 1 



Kulka- 



or ex- 



2012j) , providing evidence tor at leas t one 
UVn formation channel (see Sec 



2.11. 



Brown et al. 
plausible A1V 

In this paper, we consider five scenarios for AlVl CVn 
formation that span plausible populations, and we study 
their consequences for space-based GW detectors by 
modeling the detection and subtraction of sources in 
simulated data streams. We consider two representa- 
tive detectors: the standard LI SA developed as a joint 



NASA-ESA project until 2011 (Prince et al 



2009 Jen 



nrich 2009), and a descoped configuration with shorter 



armlengths of 1 Mkm, which is similar, but not equiva- 
lent, to the NGO/eLISA de sign studied by ESA in 2011 
( Amaro-Seoane et al.||2012 |. 

We find, first, that the number of detectable AM CVn 
systems decreases from ^ 10, 000, as reported in previous 
assessments of LISA science, to ~ 100-400 (for one year 
of observation with a single interferometer — see App.[B|). 
Precise numbers depend on the assumptions used to ex- 
plain how a postulated number of DDWD progenitors 
can result in the number of AM CVn systems seen in 
surveys. For the shorter detector, these numbers are re- 
duced by a factor of a few. Second, we confirm that 
mostly DDWD systems determine the residual confusion 
foreground left after GW signals from detectable binaries 
are subtracted. The foreground does not change signifi- 
cantly from previous estimates, except in one scenario in 
which we explain the fewer observed AM CVn systems 
by reducing the space density of DDWD progenitors. 

The rest of this paper is structured as follows. In Sec. [2] 
we discuss the current understanding of AM CVn forma- 
tion. In Sec. [H we outline the observations of AM CVn 
and DDWD systems in the Galaxy, and the resulting cal- 
ibration of local space densities. In Sec. |4]we present our 
analysis: we intr oduce our explanatory formation sce- 
narios in Sec. 



4.1 



discuss the spac e-bas ed observation of 
GWs from Galactic binaries in Sec. |4.2| describe our sim- 
ulations of GW source detection and subtraction in Sec. 
4 3| and detail our results in Sec. [5} In Sec. [6] we recapit- 
ulate our conclusions and propose possible future work. 
In the appendices we provide more detailed discussions of 



specific topics: in App. 



. X]we summarize previous analyt- 
ical and numerical studies of the Galactic foreground for 
LISA, and give a general semi-analytical derivation of its 
magnitude and shape; in App. [B] we discuss GW detec- 
tion and define our reference detector configurations; in 
App. [C] we provide more details about the fits of the GW 
confusion foreground obtained from our simulations. In 
the rest of the paper, we use the standard notation where 
G is the gravitational constant and c is the speed of light. 

2. THE FORMATION AND EVOLUTION OF COMPACT 
BINARY SYSTEMS IN THE GALAXY 

In the standard picture of binary-star evolution, com- 
pact binaries form following two consecutive mass- 



transfer episodes; the loss of sufficient angular momen- 
tum by friction in the second common-envelope phase 
ensures sufficiently tight orbits for GW emission to in- 
fluence dynamics. Here we concentrate on the processes 
leading to mass transfer and to the formation of AM 
CVn systems. Due to space limitations, we do not at- 
tempt here to prov ide an extensiv e literature review, and 
refer the reader to Solheim (20101 and references therein 
for a detailed review on the current understanding of the 
formation and evolution of DDWD and AM CVn sys- 
tems. 

The subsequent orbital evolution of compact binaries 
is driven by gravitational radiation, by mass transfer, by 
a dissipative coupling that can feed the angular momen- 
tum added through accretion back into the orbit, and by 
other poorly modeled processes such as wind loss, etc. 
Therefore, the rate of change J or b of the orbital angular 
momentum comprises three or more competing compo- 
nents, 



•/orb = Jgr. + Jm + J% 



(1) 



Here Jgr represents the loss of angular momentum to 
gravitational radiation, 



Jgr — - 



32 G 3 M X M 2 M 



orb 



(2) 



(for circular orbits) where M\ and Mi are the two binary 
masses [Mi is the donor's mass for mass-transferring sys- 
tems), a their separation, and M = M\ + Mi- Next, 
J^j = ^/GMiRhM2, where Mi indicates the donor's 
mass-loss rate due to Roche-lobe overflow and Rh de- 
fines the orbital radius around the accretor with the same 



specific angular momentu m as the transferred mass (|Ver- 
bunt fc Rappaport||l~988l ) 



Mass transfer can proceed eit her via an accretion disk 
or via "direct impact" (see, e.g., ]M arsh &; Steeghs||2002| 
Ramsay et al.|2002||Marsh et al.|2U04p . "Direct impacr 
occurs when the donor's mass stream hits the accreting 
star directly. The particular route depends on the initial 
separation in the binary following the ejection of the sec- 
ond envelope, and on the donor's entropy and equation 



of sta te (see, e.g., Gokhale et al. 2007 Dan et al. 



2011 



2012 1. Typically, mass-transferring binaries with initia 



orbital periods greater than ten minutes form an accre- 
tion disk. As discussed below, the presence or absence 
of an accretion disk is critical to the stability of the sys- 
tem. Whether a disk is present or not, J Tg denotes the 
change in the torque due to dissipative coupling, tidal or 
magnetic, between the accretor 's spin and the orbital an- 
gular momentum during the accretor's synchronization 
timescale r s . We assume that the donor and orbit's spins 
are synchronized. In this work, we do not consider the 
change in angular momentum driven by other physical 
processes, such as losses due to winds or magnetic brak- 
ing, because of current limitations in our understanding 
and modeling of such phenomena. 

As the term "detached" in their very name implies, 
DDWD dynamics are deter mined solel y by gravitational 
radiation and tides (see e.g. Piro|20lT ). By contrast, the 
formation, stability, and outcome of AM CVn systems 
depend on all three effects in Eq. (fTl). Three proposed 
channels exist for the formation of AM CVn systems, and 
each can vary in its efficiency. We describe them below. 
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Following Nelemans et al. ( 2001a|b 2004), we assume 
that the first two channels listed below can operate with 
either "optimistic" or "pessimistic" efficiencies, resulting 
in different AM CVn numbers and period distributions. 

2.1. The DDWD channel 

A subclass of DDWDs with periods of several min- 
utes becomes semi-detached when angular-momentum 
loss due to gravitational radiation drives the two bodies 
sufficiently close that the lower-mass WD fills its Roche 
lobe and mass transfer begins, either via direct impact 
or via an accretion disk. The stability of mass transfer 
depends critically on the mass ratio of the two initial 
WDs (requiring typically q = M^/Mi < 2/3) and on 
whether spin-orbit tidal coupling feeds angular momen- 
tum accreted by the primary star back to the orbit; the 
strength of this coupli ng is characterized by the synchro- 
nization timescale r s ( |Marsh et al.||2004 ). 

Under optimistic assumptions, strong tidal coupling 
stabilizes mass transfer, favoring the formation of sta- 
ble AM CVn systems and resulting in fewer mergers. By 
contrast, under pessimistic assumptions there is no tidal 
coupling and many more DDWD systems merge. In the 
presence of a disk, the accreted material's angular mo- 
mentum is mainly transferred back into the orbit by tidal 
forces at the outer edge of the disk. However, the disc's 
material at the inner edge still has sufficient angular mo- 
mentum to spin up the accreto r. Effects not considere d 
here, such as magnetic braking ( Farmer fe Roelofs|2010 ), 
may also impact the formation ot AM (JVn systems and 
favor the merger of the WDs. 

If stable mass transfer proceeds, the donor WD loses 
mass and expands in size, while the system evolves to 
longer orbital periods. The minimum orbital period at 
which mass transfer can begin is ~ two minutes for a 0.45 
M Q and 0.25 M@ binary. Maximum orbital periods for 
evolved systems in the Galaxy are estimated to be about 
90 minutes; these correspond to very old systems whose 
orbits have slowly evolved over a Hubble time. 

Surveys specifically tar geted at low-mass WD [such as 
the ELM survey of WDs: |Brown et al.| (|2010[) ; |Kilic et al. 



(2011a I; Brown et al.| ( |201l| j have resulted in the discov- 
ery of numerous DDWD systems with q < 2/3, which 
are expected to get into contact within a Hubble time. 
Several of these newly discovered low-mass binaries will 
form stable mass-transferring AM CVn systems, even in 
the most p essimistic case of no tidal coupling (Brown 
etaL|Ml|. 



2.2. The Helium star-white dwarf channel 

In this channel, a WD accretes from an initially non- 
degenerate He star with mass 0.4-0.6 Mq before the He 
has been exhausted in the He-star core. He stars are pro- 
duced from stars with masses greater than 2 Mq that 
lose their envelopes on the red-giant branch, and typi- 
cally have lifetimes comparable to their main-sequence 
progenitors, allowing for mass transfer before evolution 
to the WD stage. Mass transfer is stable when the mass 
ratio of the donor He star to the accretor WD is less than 



approximately 1.2 dTut ukov fe Fedorova||1989 
Fedorova||l990" 



see also INelemans et al. 2001a). Due to 



Ergma & 



The minimum orbital period reached by these systems is 
~ ten minutes. 

It remains in question whether 0.1 Mq or 0.3 Mq of He 
accreted from a He star will destroy a CO W D by edge 



lit de tonations (henceforth d enoted ELDs; see|Shen et al. 
2010 and references therein). |Nelemans et al.| ( [2001a[ ) dls^ 
tmguish between optimistic and pessimistic assumptions 
for this channel depending on whether double edge-lit 
detonations remove CO WD systems. 

2.3. The evolved Cataclysmic Variable channel 

In analogy to traditional CVs, a so-called evolved CV 
comprises a main-seque nce star which transfers m ass to 
a WD or neutron star ( Podsiadlowski et al.||200i" ). The 
difference is that mass transfer occurs as the donor star 
evolves off the hydrogen main sequence. However, it is 
thought that the evolved-CV channel contribute s only 
marginally (< 2%) to the AM CVn population (Nele- 
mans et al~]|2004 Roelofs et al. 2007b ) for systems with 
periods less than 1500s. Therefore, we do not consider 
this channel in this work, which focuses on the detectabil- 
ity of AM CVn systems by GW interferometers. 

3. OBSERVATIONS OF AM CVN AND DDWD SYSTEMS IN 
THE GALAXY 

We now turn our attention to recent population es- 
timates of AM CVn systems ba sed on observations. As 
mentioned, |Roelofs et al.| ( 2007b ) estimate the total space 
density of AM CVn systems using observations of a 
source population identified by emission lines in SDSS. 
Characterizing selection effects for the six observed sys- 
tems in SDSS is challenging. The orbital periods of 
these systems are generally shorter than predicted for the 
majority of AM CVn systems by population estimates, 
thereby requiring a significant extrapolation to estimate 
the space de nsity of the entire pop ulation. 

Table 1 in Roelofs et al. ( 2007b| presents the modeled 
and inferred AM CVn local space densities for differ- 
ent population-synthesis models. The inferred densities 
range from 1.2 x 10 -6 to 3.4 x 10~ 6 pc~ 3 , assuming a 
thin-disk population with a scale height of 300 pc for 
most systems. A comparison between the inferred and 
modeled densities immediately disfavors formation mod- 
els such as the optimistic DDWD channel, and requires 
an overall reduction in the estimated number of Galactic 
AM CVn syst ems by a factor larger than ten. |Roelofs| 
et al. (2007b) hypothesize that lower AM CVn densi- 



mass loss, He burning in the donor's core quenches and 
at masses ~ 0.2 Mq the star becomes semi-degenerate. 



ties arise predominantly due to the lower efficiency of 
the major formation channels. Later in this paper we 
consider an alternative possibility , namely that systems 
observed by |Roelofs et al. ( 2007b ) may come from a dis- 
tribution with scale height larger than conventional thin 
disks, with AM CVn progenitors born early in the history 
of the Galaxy, before the thin disk was fully developed. 

Although the number of known DDWD systems that 
will merge within a Hubble time has increased to 24, an 
approximate factor of five, in the last few years (four 
eclipsing systems are now known; the ELM survey has 
uncovered the eclipsing DDWD binary with the short- 
est known period of 12.7 minutes), a detailed popula- 
tion analysis of such systems has yet to be undertaken. 
More recently, on the basis of the assumed Galactic-disk 
model, the ELM survey predicts 40 kpc -3 of DDWDs 
where at least one WD has a mass < 0.25 Mq, with a fac- 
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tor of two uncertainty in correction factors (Brown et al. 
2011 1 . This is much lower th an population-synthesis re- 
sults ( Nelemans et al.|200Tc |. On the other hand, popu- 
lation syTrnTesIsresulTssnow good agreement with obser- 
vations made by the Supernova type la Progenitor (SPY) 
survey. The SPY radial-velocity survey estimate d the lo- 
cal space de nsity of WDs to be 4.8 x 10~ 3 pc~ 3 (Holberg 
et al.||2008|) , with ~ 100 new DDWDs found (with or- 
bital periods less than ten days) from a sample of 1000. 
At present, estimating space densities from the ELM 
and SPY surveys is very challenging, requiring proper 
understanding of selection biases. Detailed population 
comparisons between the two surveys are not currently 
available. 

4. OUR ANALYSIS 



We now describe our analysis: in Sec. 4.1 we expl ain 
our choices for WD binary populations; in Sec 
discuss GW emission from binaries; and in Sec 



4.2 



4.3 



we 
we 

describe our simulations of binary detection and subtrac- 
tion in detector data. Sec.j5] below presents the resulting 
estimates of the number ofdetectable systems and of the 
residual confusion-noise foreground. 

4.1. WD binary catalogs 

In our analysis we considered several different assump- 
tions on AM CVn formation, generated simulated bi- 
nary catalogs and detector datasets for each, and es- 
timated the resulting binary detection and subtraction 
prospects. Among the large number of cases that we 
investigated, here for conciseness we report on five sce- 
narios that we believe are representative of the range of 
possibilities given the constraints of current observations. 
Two of the scenarios represent an upper bound (with 
maximally efficient formation channels) and conversely, 
a bound derived assuming inefficient formation channels 
on predictions from population-synthesis models (cases 1 
and 2); two others attempt to reproduce SDSS-estimated 
AM CVn densities (cases 3 and 4 below); a fifth (case 5) 
explores a Galactic population with a larger scale height 
at times earlier than 7 Gyr ago. As discussed below, the 
modified-Galactic-disk scenario (case 5) should be con- 
sistent with SDSS observations; however, this agreement 
is difficu lt to quantify as the S PSS-estimated densities 



given in Roelofs et al. (2007b) are derived assuming a 



standard disk model. 

We generate DDWD and AM CVn catalogs with the 
SeBa popu l ation-synthesis code (|Portegies Zwart fe Ver- 
bmrfc] [T9Ml |Nelemans et~aLl|2001b[ \Ml\ |Toonen et, al. 
20111, with the following assumptions and parameters: 
i) the time- and position-dependent star form atio n rate 



(SFR) based onlBoissier fe Prantzos| ( |1999[ ) (see [Nele- 
mans et al.||2004 j); ii) initial pr imary masses distr ibuted 
with power-law index —2.5 (see Kroupa et al.|T9 93); iii) a 
flat initial mass-ratio distribution; iv) a flat distribution 
for the logarithm of the semimajor axis, up to a radius 
of to 6 /?©; v) eccentricities distributed as P(e) oc e 2 ; vi) 
a 50% binary fraction in the initial population of main- 
sequence stars. We expect a factor of two difference in 
population-synthesis results obtained using different as- 
sumptions and initial parameters than those stated above 



(Portegies Zwart & Verbunt 1996; Ne lemans et al.|2001b 



2004 ) . Each system is described by seven parameters 
the binary masses Mi and M 2 , the orbital period P Q rb 5 



the mass-transfer rate M, the Galactic latitude / and 
longitude 6, and the luminosity distance d to the Solar 
System. 



Following Nelemans et al. (20041, we distribute single 
and binary stars according to a Galactic model that com- 
prises both disk and bulge components]^ The density of 
systems in the disk is 



Pdisk(-R, Z) = pbp(R) sech 2 



Z 



J disk 



pc 



(3) 



where R is the cylindrical radius from the Galactic cen- 
ter; pBp(R) results fro m integrating the time - depen dent 
plane-projected SF R of Boissier fc P ran tzos ( 1999 ) [see 
Eqs. (1) and (2) of |Nelemans et al.| ( |2004| j; Z isTne per- 
pendicular distance from the Galactic plane; and Zdisk 
is the disk's scale height. The density of systems in the 
spherical bulge is given by 



Pbuigo (r) oc e 



-(''/''bulge)' 



pc 



(4) 



where r is the spherical distance from the Galactic center; 
7"buige is the characteristic radius of the bulge; and Eq. 
Q is normalized so that there are as many bulge systems 
as disk systems in the inner 3 kpc of the Galaxy, with a 
current total mass of 2.6 x 10 6 M . 
For the first f our scenarios we set i?disk = 300 pc and 
0.5kpc (Juric et aT]|2008 Nelemans i et al 



^bulec 



2004 



Binney fc Tremame|fl987| |sbtue et al.||2009| ); pbp(R) is 



ht well by a density oc exp(— /i/f^isk) with i?disk ~ 2.5 
kpc for DDWD systems and i? disk « 2.2-2.4 kpc for AM 
CVn systems. In the last scenario we explore a modified 
model of the Galactic disk with a larger scale height for 
binaries that are produced at earlier times ( > 7 Gyr ago) . 

For all bu t one scenarios we use the DDWD catalog- 
analyzed in Nelemans et al. (2004), with local density 
1.2 x 10 _5 pc _,:i , and orbital periods as large as 5.5 hrs 
(/gw = 10~ 4 Hz). Binaries at longer periods are not 
typically detectable by GW instruments because detector 
noise is higher at low frequencies; they are therefore not 
considered further in this analysis. All the catalogs of 
AM CVn systems, by contrast, show the effects of our 
different assumptions on AM CVn formation: 

Case 1: upper bound (strong tidal coupling and ineffi- 
cient ELD events). This catalog includes AM CVn 
systems that have formed via the DDWD channel 
operating with strong tidal coupling and t he He 
star- WD channel with fe w ELD events (see Nele 



mans et al. 2001b 2004 1. Although the resulting 



AM CVn local density is higher (2.8 x 10~ 5 pc~ 3 ) 
than SDSS-estimated values, past investigations of 
Galactic WD binaries as LISA sourc es have all used 
variants of this scenario (see, e.g., |Edlund et al . 

Thus, this case pro 



2005 



[Timpano et al. 2006) 
vides a meaningful comparison between our study 
and previously published results, which may vary 
depending on their ideali zatio ns of "detection," 



"resolution," etc. (see Sec. 4.3) 



Case 2: inefficient formation channels (weak tidal cou- 
pling and efficient ELD events). This catalog in- 



6 The bulge is not included in the a nalysis of iRoelofsetal. 
1 2007b) . Furthermore, both |Ruiter et al.| l [2009| ) and |Yu fc Jettery 
1 20iu[l a lso include a halo component; however, the impact ot halo 
BTTWT) systems seems negligible for GW interferometers. 
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eludes AM CVn systems that have formed via the 
DDWD channel operating with minimal tidal cou- 
pling and t he He star-WD channel with m any ELD 
events (see |Nelemans et al.|2~001b[|2004[ ). The AM 
CVn local density is 6.1 x 10 _t> pc -3 , a factor of 



two greater than SDSS-estimated values (Roelofs 



et al. 2007a I , and fall within the factor of two un 



certainty ot population-synthesis results. 

Case 3: mild tidal coupling and no ELD events. This 
catalog includes AM CVn systems that have 
formed via only the DDWD channel, where mild 
tidal coupling (r s = 2 years) between the WDs 
allows for stable mass transfer in up to 10% of 
systems compared to the optimistic channel of 
Case 1. We assume that the He star-WD chan- 
nel is suppressed. [Howev er, observations discussed 
in ( |Roelofs et al. 2007a) suggest the presence of 
hot donors, which imply either that the He star- 
WD channel is active with core He burning hav- 
ing been quenched only recently, or that follow- 
ing the second common-envelope ejection, the ini- 
tial DDWD orbital separation is so small that the 
donor WD cannot cool within the gravitational- 
radiation timescale.] The AM CVn local density is 
1.1 x 10" 6 pc~ 3 . 

Case 4: fewer progenitor DDWDs. In this scenario we 
explain the SDSS-estimated AM CVn densities by 
postulating a 1/5 reduction in the number of their 
DDWD progenitors. This reduction in DDWD sys- 
tems is consistent with existing observations and 
with the factor of two uncertainty introduced in 
population-synthesis models by our limited under- 
standing of binary evolution and Galactic struc- 
ture. This catalog also includes AM CVn systems 
that have formed both via the DDWD channel with 
mild tidal coupling (r 8 = 1 year) and via the pes 



simistic He star-WD channel (see Nelemans et al 



200Tbl 



2001b| 2004); their densities are also reduced by 



1/5 (across all periods). The DDWD and AM CVn 
local densities are 2.4 x 10~ 6 and 1.5 x 10~ 6 pc~ 3 . 

Case 5: modified Galactic-disk model. In this scenario 
we place older systems in a thicker disk. With the 
caveat of small-number statistics, AM CVn sys- 
tems identified by emission lines in SDSS have in- 
ferred projected velocities > 70 km s , hinting 
that AM CVn systems with longer periods may 
originate from a disk with a thick component, or 
at least from a thin disk with a larger scale height. 
Indeed, such a distribution of AM CVns at birth 
is consistent with our current (if everchanging) 



tion of thin and thick g 


alaxy disks flScannapieco 


et al. 120111 
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ited (see |Napiwotzki||2009| for a discussion on the 
thin-disk, thick-disk, and halo populations of single 
WDs). As an illustrative example, we altered the 
disk component of Eq. [3] to have Zdisk = 1250 pc 
for systems older than 7 Gyr, Z^isk = 300 pc for 
systems younger than 7 Gyr, and rtmige = 0.5 kpc 



as for cases 1-4. This effectively places all sys- 
tems older than 7 Gyr in a thick disk, which in- 
cludes the great majority of longer-period systems. 
In addition, we assumed that AM CVn systems 
only formed via the DDWD channel. The DDWD 
and AM CVn local densities are 1.2 x 10~ 5 and 
3.8 x 10~ 6 pc- 3 . 

For each scenario, Table[l]shows the numbers of AM CVn 
and DDWD systems in the catalogs, as well as their local 
density po (i.e., the Galactic mid-plane density at 8.5 kpc 
from the Galactic center). 

4.2. Detecting gravitational waves from WD binaries 

Compact Galactic binaries in the low-frequency band 
emit quasi-monochromatic GWs, with frequency drifts 
of (at most) several frequency bins for year-long obser- 
vation. In this paper, we model compact-binary wave- 
form s by including only the d ominant quadrupolar emis- 
sion ( Peters fc Mathews|1963 ), with instantaneous strain 
amplitude 



lO" 21 



Mz_ 



5/3 / p N -2/3 
r orb 

hr 



d 
kpc 



(5) 



where M z = (1 + z)/i 3 / 5 M 2 / 5 = (1 + z)M is the red- 
shifted chirp mass (with M the total mass, the reduced 
mass, M. the chirp mass, and z the redshift), d is the lu- 
minosity distance to the source, and P rb is the orbital 
period of the binary. 

Figure [T] shows the number density as a function of h 
[the GW strain of Eq. @], and / = 2/P orb [the GW 
frequency], for the different DDWD and AM CVn cata- 
logs of each of our scenario s. T he colored dots indicate 
detected binaries (see Sec. 4.3): green shows DDWDs, 



red shows AM CVn systems (from the DDWD channel) 
and blue shows AM CVn systems (from the He-star- WD 
channel) . 

Binaries with orbital periods shorter than approxi- 
mately 20 minutes will appear as isolated signals in year- 
long datasets of space-based GW observatories, and will 
be detectable when the signals are sufficiently strong 
compared to instrument noise at the same frequency. 
By contrast, most binaries with longer orbital periods 
will not be detected individually, because they will be 
too dense in frequency space to tell apart. These will 
produce a noise-like confusion foreground that will affect 
the detection and parameter estimation of extra-Galactic 
GW sources such as massive-black-hole binaries. Never- 
theless, the brightest binaries will be detected above the 
confusion foreground. 

Looking at Fig. [l] we see that DDWD systems have 
larger densities than AM CVn systems, and larger GW 
strains (because they have greater chirp masses). Thus 
it is DDWD systems that are principally responsible 
for the shape and stre ngth of the GW foreground from 



WD binaries (see Sec. 5.2). We note also that all AM 
CVn populations show an overall low-frequency cutoff at 
log 10 / = —3.5, due to the age of the Universe. 

A simplified (but roughly correct) description of GW 
data analysis is that GW signals will be detectable when 
their signal-to-noise ratio, SNR, is greater than a detec- 
tion threshold SNR t hr chosen to yield a minimal number 
of false alarms. SNR is the ratio of signal strength to 
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TABLE 1 

Total number and local density of compact-binary systems of different classes (detached, AM CVn from the DDWD 
channel, AM CVn from the He-star-WD channel) in our simulated catalogs for cases 1-5. 

number in catalog / local space density [pc — 3 ] 
DDWD AM CVn (DDWD ch.) AM CVn (He-star-WD ch.) 



case 1 
case 2 
case 3 
case 4 
case 5 



26,084,411 / 1.2 X 1CT 5 
26,084,411 / 1.2 x 10 -5 
26,084,411 / 1.2 x 10~ 5 
5,217,866 / 2.4 x 10~ 6 
26,084,509 / 1.2 x 10~ 5 



23,025,764 / 1.9 x 10~ 5 
261,840 / 2.1 x 10~ 7 
3, 178,553 / 1.1 x 10~ 6 
994, 194 / 2.6 x 10~ 7 

23,025,764 / 3.8 x IP' 6 



11, 197,735 / 9.0 x 10 
6,643,091 / 5.9 x 10 

none / 
1,328,909 / 1.2 x 10 

none / 



case 1 



case 2 



case 3 



case 4 



case 5 



d/V 



dlog 10 /7dlog 10 f H; 
DDWD 



AM CVn 

(DDWD and He ch.; 
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Fig. 1. — Number density as a function of GW-strain amplitude h and GW frequency / for DDWD (top row) and AM CVn (bottom 
row) systems in our five scenarios. Green dots denote individually detected DDWD systems (for one year of observation with the 5- Mkm 
detector); red dots AM CVn systems from the DDWD channel, and blue dots AM CVn systems from the He star channel. See Sees. |473] 
and |5.1| for a discussion of individually resolvable sources. The noticeable low-frequency cutoff of AM CVn systems is due to the age of the 
Universe. 



rms instrument noise (plus confusion noise, if present) 
at the same frequency. For LISA-like interferometric 
detectors, the yearly orbit and rotation of the space- 
craft constellation imprint the signals with frequency 
and amplitude modulations that encode the sky position 
of the source. Detector response is further affected in 
frequency-dependent fashion by Time-Delay Interferom- 
etry (TDI) , the technique used to subtract the otherwise 
overwhelming laser frequency noise (thus, SNR would be 
evaluated for the TDI signal with respect to TDI noise) . 
See App. [Blfor an overview of LISA-like detectors, their 
response to GWs, and the theory of probabilistic signal 
detection in noise. 

4.3. Simulations of detection and subtraction 

To estimate the number of Galactic-binary systems 
that we will be able to detect individually with our ref- 
erence detectors, as well as the confusion noise from 
the residual unsubtracted systems, we performed sim- 
ulations of detection and s ubtraction using a varia nt of 
the met hod introduced by Timpano et al. (2006) (see 



2007) 
a tul. 



Sec. 5.1). To wit: 



1 



For each of our scenarios we compute the 
frequency-domain TDI response for each system 
in the c orresponding ca talog, using an optimized 
variant ( |Vallisneri||2011[ ) of the Mock LISA Data 
Challenge fast- binary code (Cornish & Littenberg 



; we then superpose the responses to produce 
year-long dataset. 



2. We run through the catalog and evaluate the SNR 
of each system (see Sec. [B]) with respect to a back- 
ground of instrumental noise plus confusion noise 
from all nearby systems, estimated as a running 
median of their power spectral density. We flag as 
detected the systems with SNR > SNR thr . 

3. We subtract the signals detected in step 2 from the 
dataset. We then go back to step 2 and iterate, us- 
ing the resulting lower estimate of confusion noise. 

4. The process converges after a few iterations, with 
the number of detected systems decreasing better 
than exponentially. We then evaluate and fit the 
power spectral density from the residual confusion 
foreground. 

This process represents a rather optimistic model of 
source detection and subtraction that neglects two im- 
portant effects: source detection will be confused by the 
presence of nearby signals beyond what is predicted by 
the increased noise level, and subtraction will be de- 
graded by the imperfect estimation of parameters in noise 



(see, e.g., Crowder & Cornish 2004 1. Nevertheless, we 
can adjust the level of optimism by varying our assump- 
tions regarding the duration of observation T Q b s , the de- 
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tection threshold SNRthr, and the number iV b s of TDI 
observables that are available (two, in effect, for the stan- 
dard LISA configuration with five or six operating links; 
one for a two-arm configuration) . 

It turns out that our results for the number of detected 
sources and for the residual confusion foregrounds are fit 
remarkably well by power laws of a single effective-SNR 
parameter p cS = SNRth^Tobs/yr)- 1 / 2 ^^)- 1 / 2 . As a 
baseline we take a standard assumption, consistent with 
prior literature, of one year, one observable (the standard 
unequal-arm Michelson X), and SNRthr = 5. A conser- 
vative assumption would correspond to taking one year, 
one observable, and SNR th = 10, yielding p cS — 10; an 
optimistic assumption to five years, two independent ob- 
servables, and SNRth = 5, yielding p c s = 1.58. Results 
are not sensitive to details such as the choice of running- 
median window, and whether signals are subtracted im- 
mediately or after each full catalog run-through. 

5. RESULTS 

In this section we present our results on the impact that 
different Galactic populations of DDWD and AM CVn 
systems have on the number of individually detectable 
sources (Sec. |5.1|) and on the total and confusion GW 



foregrounds (Sec. 5.2), for our two representative space- 



based GW detectors and for each of our five population 
scenarios. 

Figure [2] recapitulates the findings of our simulations: 
it shows the frequency-space density of catalog sources; 
the power spectral density of the unsubtracted (blue) and 
subtracted (red) confusion foreground, compared to the 
instrument noise (black, all in the bottom panel); and 
the location and numbers of detected sources (dots in 
the bottom panel, histograms in the top panels). The de- 
tected sources lie approximately a factor of SNR 2 h above 
the composite of instrument noise and subtracted confu- 
sion noise. 

5.1. Individually detected sources 

Tables [2] (for the 5 Mkm detector) and [3] (for the 
1 Mkm detector) show the number of recovered sys- 
tems of each type in our five scenarios. The num- 
bers quoted correspond to simulation results for one 
year, one interferometric observable, and SNRthr = 5; 
the power laws give approximate fits as functions of 

Pes = SNR th r(Wyrr 1/2 (AU s r 1/2 . 

Consider first the number of detectable sources as seen 
by the longer, standard-LISA-like detector. In all cases 
except the "fewer progenitors" scenario 4, the number 
of detected DDWDs remains comparable to prior opti- 
mistic projections (~ 12,000-14,000). Even in case 4, 
where the DDWD systems in the catalog are reduced by 
a factor of five, detections drop only by half, thanks to 
a weaker confusion foreground. By contrast, the number 
of detected AM CVn varies considerably depending on 
their density and formation scenarios (~ 15-2,000). In 
cases 2 and 3, more DDWDs are detectable at relatively 
high frequencies because of the reduced confusion noise 
from fewer AM CVn. 

Case 1 was chosen to be directly compara ble to the 



standard population ca talogs of earlier works ( Nelemans 

3er of de- 



et al. 2001a b c 



tections for A 



2004 1, but the claimed num 
CVn systems is significantly lower than 



those earlier estimates, derived with a simple resolvabil- 
ity criterion of a single so urce per frequency bin with 
instrument-noise SNR > 1 ( |Nelemans et al.|200TB]|2004[ ). 
This is unsurprising, given that in our analysis we model 
detection by identifying and subtracting bright sources 
iteratively, and that we require a higher SNR, defined 
with respect to both instrument and confusion noise. 
Furthermore, in this paper we use proper interferometric 
observables instead of "raw" GW strain. 

Downscaling the detector's armlength to 1 Mkm ap- 
proximately halves the number of DDWD systems de- 
tected for p c ff = 5, except for case 4, where the decrease 
is ~ 3. The number of AM CVn system detections is 
likewise reduced by factors ~ 3-4. Note however that 
the p e fi scalings of Tables [2] and j3] are slightly different . 
The ESA eLISA/NGO study ( |Amaro-Seoane et al |2012 1 
finds numbers compatible with ours, once we take into 
account our noise assumptions for the 1 Mkm detector, 
which are somewhat more favorable than NGO's. 



5.2. Unsubtracted and residual compact-binary 
foregrounds 

In Appendix |Aj we apply a familiar analytical treat- 
ment of the compact-binary GW foreground to a broad 
range of Galactic models examined in the literature, and 
we conclude that the foreground depends strongly on 
the Galactic binary merger rate and on their character- 
istic chirp masses, but only weakly on other parameters. 
Namely, 



S h (f) ~ 1.9 x 10- 44 (//Hz)- 7 / 3 Hz- 1 



(6) 



V, 



char 



11 



gal \ 



6.4 kpc/ V°- 015 /y i 7 V°- 35M < 



M 



z, char 



5/3 



where 2? c har is the characteristic distance from Earth to 

Galactic binaries (a weak function of the model); 

is the binary merger rate in the Galaxy; and A4 ZtC h a r 

is the characteristic chirp mass (A4^ 3 )~ 5 / 3 . The /~ 7 / 3 
slope describes the power emitted at GW frequency / 
by a single binary evolving due to GWs, and it applies 
to a stationary ensemble where new binaries are born 
with constant rate (equal to 7?. ga i) at frequencies below 
the detector bandwidth. Equation ^ is normalized to 
values representative of our assumptions; as shown in 
the left panel of Fig. [3j it is in good agreement with the 
results of our simulations for cases 1-3 and 5, while the 
coefficient is reduced by a factor of five for case 4. 

The residual confusion foreground consists of the un- 
resolvable DDWD and AM CVn systems that are left 
once relatively brighter or isolated sources have been sub- 
tracted iteratively from the original population catalog. 
The residual foregrounds found in our simulations are 
very similar among cases 1-3 and 5, while they are a 
factor of five weaker for case 4. We note that the over- 
all number of DDWD systems in cases 1-3 and 5 differ 
by a factor of five, thereby confirming earlier work and 
our predictions that the DDWD population provides the 
dominant component to the strength of the confusion 
foregrounds (both unsubtracted and subtracted). As 
shown in the middle and right panels of Fig. [3j for both 
the 5 Mkm and 1 Mkm detectors, the total equivalent 
strain noise (instrument noise plus confusion foreground) 
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Fig. 2. — Frequency-space density and GW foreground of DDWD and AM CVn systems for the five scenarios examined in this paper. 
In each subplot, the bottom panel shows the power spectral density of the unsubtracted (blue) and partially subtracted (red) foreground, 
compared to the instrument n oise (black). The open white circles indicate the frequency and amplitude of the so-called "verification 
binaries" (see [Nelemans 2012). These comprise known DDWD and AM CVn systems; a handful have known distances from HST FGS 
parallax observations, while we assume that the rest are at a distance of 1 Kpc. The green dots show the individually detectable DDWD 
systems; the red dots show the detectable AM CVn systems formed from DDWD progenitors and the blue dots show the detectable AM 
CVn systems that arise from the He-star- WD channel. The top panel shows histograms of the detected sources in frequency space. 
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TABLE 2 

Detectable sources in our simulations (5 Mkm detector). Values are shown for SNR thr = 5, one interferometric 

OBSERVABLE, AND ONE YEAR OF OBSERVATION, WITH APPROXIMATE SCALINGS AS A FUNCTION OF p e ff = SNR t hr(7 1 b s /yr) _1 / 2 (A r bs)~ 1 ' /2 - 







DDWD 


AM CVn (DDWD) 


AM CVn (Hc-WD) 


case 


1 


11,898 x (peff/5)- 1 - 4 


2,093 x (peff/5)- 1 - 8 


122 x (p eff /5)- 2 - 8 


case 


2 


13,957 x (pcft/5)- 1 - 4 


16 X (Peff/5)" 1 - 7 


72 x (p cff /5)- 2 - 9 


case 


3 


13, 708 x (peff/5)- 1 ' 4 


397 x (p off /5)- 2 - 


none 


case 


1 


5,252 x (peff/5)" 1 - 2 


161 X (Peff/5)" 1 - 7 


79 X (Peff/5)" 1 ' 8 


case 


5 


11,972 x (pea/5)- 1 - 4 


2,107x (pes/5)- 1 - 7 


none 



TABLE 3 

Detectable sources in our simulations (1 Mkm detector). Values are shown for SNR thr = 5, one interferometric 

OBSERVABLE, AND ONE YEAR OF OBSERVATION, WITH APPROXIMATE SCALINGS AS A FUNCTION OF p ef f = SNR thr (T obs /yr)~ 1 / 2 (A r obs ) -1 / 2 . 







DDWD 


AM CVn (DDWD) 


AM CVn (Hc-WD) 


case 


1 


6,034 x (peff/5)- 1 ' 2 


719 X (peff/5)- 1 ' 4 


38 x (p eff /5)- 2 - 6 


case 


2 


6,337 x (poff/5)- 1 ' 2 


4 X (Peff/5)" 1 - 4 


25 x (p cff /5)- 2 - 2 


case 


3 


6,306 x (peff/5)- 1 - 2 


87 X (peff/5)" 1 ' 8 


none 


case 


4 


1,551 x (peff/5)- 1 - 1 


40 X (peff/5)- 1 ' 3 


8 x (p cff /5)- 2 - 2 


case 


5 


6,155 x (poff/5)- 1 - 2 


727 x (pes/5)- 1 ' 4 


none 




f[Hz] f[Hz] f[Hz] 

Fig. 3. — Left panel: power spectral density of GW strain for the unsubtracted Galactic foreground, as measured in our simulations 
for cases 1—3 and 5 (top continuous curve) and case 4 (bottom continuous curve); the foreground is fit well by our Eq. I |6|l, w ith a 1/5 
smaller coefficient for case 4 (dashed curves). For comparison, the dotted curves shows the Barack and Cutler estimate Eq. |[C2l. At high 
frequencies, the simulated spectra drop because frequency bins are unevenly occupied, and because of the onset of mass transfer; at low 
frequency, they drop because some binaries are born at frequencies within the detector's bandwidth, contrary to the assumptions that yield 
the f~ 7 / 3 profile. [These spectra of GW strain do not include the instr umen t response to GWs; to compare with the TDI X spectra of 
Fig. [2] one would multiply by the frequency-dependent response of Eq. 1C8L ] Middle and right panels: total equivalent GW-strain 
noise(instrument noise plus residual Galactic-binary confusion noise) for the 5-Mkm and 1-Mkm detectors of our study, as obtained in 
our simulations for cases 1—3 and 5 and for p^g = 5 (red curv es); c ase 4 is not shown to avoid visual clutter. The total noise is fit well 
(dashed curves) by our expressions Eqs. Q- (j9Jl [see also Eqs. ( |C4[ )— |C7| ] , and less well (dotted curves) by Barack and Cutler's standard 
prescription Eqs. l|Cl|)— (|C3[I, which improves by setting kT" 1 ~ 9/yr. 
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is fit well by the expression 



rrtotal 



(/) = a*-* + S c h oni x tanh° 



[tfdN/df), (7) 



where a ~ 1, /3 ~ 1/yr, 



^° nf (/) = 1-4 x 10- 45 (//Hz) 8/3 Hz" 1 , (8) 
and where 



dV 



= 5x 10- 3 (//Hz) _11/3 Hz- 1 



(9) 



is the frequency-domain density of binary systems, which 
modulates the transition from the confused to the re- 
solved regime. Equation J7J is fo rmula ted in analogy to 
the stan dard expression Eq. (CI) [see (Barack & Cutler 
|2004a|bj )] , with a different functional form that tits our 
simulations better. In App.[C]we give more precise a and 
/3 fits as functions of p e ff > an d we provide expressions for 
the total TDI X noise and for the seasonal mo dulation of 
the GW foreground due to detector motion dG iampieri 



fc Polnarev|p97l [Setol[2004l |Edlund et al1[200^ 



6. CONCLUSIONS 

In this work, we have examined the impact that dif- 
ferent observationally inferred populations of Galactic 
white dwarf binaries have on the number of individu- 
ally detectable sources and on the confusion foreground 
for future low-frequency GW observations. Earlier stud- 
ies of Galactic white dwarf binaries as GW sources re- 
lied primarily on population synthesis, because of the 
small number of observed systems, in particular mass 
transferring ones (i.e., AM CVn systems). However, the 
past several years have seen a flurry of new observations, 
quintupling the known detached systems that will merge 
within a Hubble time, doubling the known AM CVn sys- 
tems, and leading to new, lower estimates for the local 
densities of the latter. 

To understand the consequences of these new estimates 
on GW detection, we have generated a range of plausi- 
ble populations of white dwarf binaries, assuming differ- 
ent AM CVn formation scenarios. For example, in one 
model (case 3) we postulate that fewer AM CVn sys- 
tems are observed by SDSS than predicted because the 
tidal coupling between the spin of the accretor and the 
orbit is not as efficient as expected, so more systems are 
lost to mergers. We have then investigated the resulting 
numbers of individually resolvable GW sources and the 
magnitude of the total and residual GW foregrounds. 

In four out of five models, we find that the number of 
detectable detached systems is ~ 12,000-14,000, for one 
year of observation by a LISA-like GW detector with 5 
Mkm arms (using an optimistic detection threshold, but 
considering a single interferometric observable). These 
four models leave unchanged the number of detached bi- 
naries used in earlier estimates, but make various as- 
sumptions about the number of systems lost to merg- 
ers or to edge-lit detonations, or (in case 5) about the 
spatial distribution of binaries. In one model (case 4), 
which instead explains the observed local AM CVn den- 
sity with five times fewer detached progenitor systems, 
we find that ~ 6,000 of the latter should be detectable. 
For the 1-Mkm-arm detector, DDWD detections would 
be ~ 6,000 and 1,500 for cases 1-3 and 5, and case 4, 
respectively. The number of detectable AM CVn systems 



ranges from several tens to a few hundred (for cases 2-4), 
or even as high ~ 2,000 (for cases 1 and 5), again for the 
5-Mkm detector. For the 1-Mkm detector these numbers 
are smaller by approximately a factor of three to four. 
See Tables |2] and [3] for details. 

Now, why is it that even in our upper-bound model 
(case 1), where the total number of AM CVn systems 
is comparable to earlier estimates, we projec t far fewer 
AM CVn detection s than the 10,000 quoted by Nelemans 
et al.| ( |2001b"l|2004| )? Conversely, in four out of five mod- 
els we project more det ached-binary detections ( 12,000- 
14,000 vs. 10,000) than |Nelemans et al.j ( [200Tb| . Such 
differences arise because we define detection by way of an 
iterative identify-and-subtract process, where a source's 
signal-to-noise ratio is evaluated with respect to noise 
from both the instrument and the partially subtracted 
foreground. Our app roach is similar to that us ed by 



Imans et al. "(2001b 



Timpano et al. (2006). By contrast, according to Nele 



]2004[ ), a source is detectable if it 
has a signal-to-noise ratio (relative to instrument noise 
alone) greater than five, or a signal-to-noise ratio greater 
than one with no other source in the same frequency 
bin. Our more stringent definition impacts the number 
of detected AM CVn systems more than their detached 
counterparts at the same frequencies, because AM CVn 
systems have smaller chirp masses, and therefore rela- 
tively weaker GW-strain amplitudes. 

Indeed, we confirm that the dominant contribution to 
the diffuse GW foreground comes from detached bina- 
ries, which have larger chirp masses. In App. [XJ we 
show that the spectrum of the unsubtracted foreground 
can be derived robustly by simple analytical arguments, 
and depends mainly on the Galactic merger rate and the 
characteristic binary chirp mass [see Eq. ([6])]. These an- 
alytical estimates agree well with our numerical results. 
In addition, for both the 5 Mkm and 1 Mkm detectors we 
provide fitting expressions [Eqs. ([7])-([9|] for the residual 
confusion foreground due to unresolvable sources. Our 
estimates are comparable to (but meas urably different 
fro m) the s tandard formulas discussed by Barack & Cut- 
ler] ( 2004b ) and used by many authors in the past. 

The projected detection numbers and residual con- 
fusion foregrounds change significantly if we use dif- 
ferent assumptions on the detection threshold, the du- 
ration of observation, and the number of available 
interferometric observables. We have simulated the 
identify-and-subtract process for several such assump- 
tions, and we find empirically that the detection num- 
bers and confusion-foreground fitting parameters scale 
as power laws of a single effective threshold p c g = 
SNR thr (r obs /yr)- 1 /2( 7Vobs )-i/2_ 

In one of our models (case 5), we have examined (to 
our knowledge, for the first time) how GW observations 
would change if we modify the evolution and structure of 
the Galactic disk to match the observed local AM CVn 
density, while keeping overall binary numbers as in case 
1. We find no significant impact on the number of de- 
tectable systems, nor on the strength or shape of the 
confusion foreground. However, on average the detected 
systems will have measurably smaller GW amplitudes, 
because of their larger distances from the Solar system. 
We conjecture that such distance measurements could be 
used to infer the structure of the Galactic disk. 
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TABLE 4 

Estimates of GW from Populations of Compact Binaries 



reference 



model 



source classes 



input population 



(2004) RS + CB DDWD + AMCVn ScBa 1 



this work, Nclemans et al. 



Nelemans et al.| ( 2001b ) 
Edlund et ail, j2005| 



RS 
RS 



CB 
CB 



DDWD 
DDWD 



AMCVn 
AMCVn 



Timpano et al. ( 2006 ) 
Seto 1 2002 ) 
B i & Hils (1997) 
Hils et al.|fn$90t 



RS + CB 
RS + CB 
RS + CB 
CB 



DDWD 
DDWD 
DDWD 
DDWD 



ScBa 1 

SeBa 1 

ScBa 1 , Wcbbink 1 1984) 

Wcbbink (1984) ~_ 

~ Wcbbink 1 1984) 
Webbink (1984) 



Hils & Bender ( 2000 ) 
Hils||l998|) 



CB 
RS + CB 



AMCVn 
AMCVn 



Ibcn & Tutukov ( 1984 ) 



Yu fc Jcff cry ( 2010 ) 

Liu|||20 09); Liu et al. (20W)_ 



Tutukov fc Yungelson 



iben & 'l\itukov||i984l; 'Ibtukov fc Yungelson 



RS + CB DDWD (+ halo) 



Farmer & Phinncy ( 2003 1 



RS + CB DDWD 



BSE 2 
BSE 2 



iHanj (1998) 



(1996) 



(1996) 



CB extragalactic DDWD BSE 2 



- Han ( 1998 1 



+ AM CVn 


Wcbbink & Han (1998 


) CB DDWD 


Han 


(1998 


) 



Ruiter et al.| | |2010[ | 
Ruiter et afl J2OO9I 



RS + CB 



DDWD + AMCVn 
(halo) 



StarTrack, Bclczynski ct al. ( 2002 2008 1 



Benacquista et al. (2004) 



RB + CB DDWD (+ globular) Specified in paper 



Postnov fc P rokhorov I1998J 
Lipunov et al. ( 1995 ) 
Lipunov & Postnov (1987) 

Lipunov et a l. (1987) 

Giampicri & Polnarcv (1997) 



CB 
CB 
CB 



DDWD 

DDWD (+ halo) 
compact binaries 



See |Lipunov et al.| | |1996| 

specified in paper 
specified in paper 



Evans et al.l 111987 1 



CB 



generic 



gene 1 lc 



CB 



DDWD 



specified in paper 



Notes: RS = resolved sources; CB = conf u sion background. 
1 SeBa : [ Portegies Zwart fc Verbu nt 
2 BSE: iHurley et all \2000\ |2002T" 



119961; Portegies Zwart & Yungelson 1 19981; Nelemans et al. ( 2001c 1; Toonen et al. 1 2011 1 



Finally, we expect that the advent of wide-field synop- 
tic and targeted surveys (such as the Palomar Transient 
Factory and the ELM survey respectively), as well as 
ESA's GAIA mission, will continue to provide new data 
to this field of research, and will enable more stringent 
constraints on binary populations. So far, the observed 
density of detached systems appears to be consistent with 
theoretical estimates, but a detailed population analy- 
sis based on the ELM survey has yet to be undertaken. 
Furthermore, all estimates of AM CVn densities based 
on surveys are made uncertain by the intricate selection 
effects. Because of these uncertainties, in this paper we 
chose to examine a variety of population models that 
span the range of likely outcomes. Our overall conclu- 
sion is that the Galactic-binary science enabled by low- 
frequency GW detectors is very robust, for a broad range 



of possible binary populations and detector sensitivities. 
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APPENDIX 

PREVIOUS ESTIMATES OF GW EMISSION FROM GALACTIC BINARIES 

More than two dozen papers have made estimates of the GW signal from individually detectable Galactic compact 
binaries and from the continuum of unresolvable sources at low frequencies (variously called "confusion background" or 
"Galactic f oreground" ) for a LI SA-like GW detector with 5-Mkm arms. An earlier informative review of the literature 
is given by Ruiter et al. (2010). We note that for a 1- Mkm detector the foreground decreases substantially, to a level 
comparable to instrument noise (see Fig. [3] in Section 5.2). 

Table H compares a subset (large, but possibly not complete) of relevant published estimates, grouped by population 
model. For each entry, we show whether the estimates include resolved sources (RS), the continuum background (CB), 
or both; we show also which classes of binaries are included. In Table [5] we provide relevant scale lengths for the 
Galactic disk and bulge distributions, where available; these are indicative of the range of assumed Galactic models, 
but they do not provide a complete picture since different density laws also appear. We do not give quantitative 
information on possible halo components, which are also frequently used. The major differences between the various 
estimates arise from two sets of assumptions: 
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TABLE 5 

Estimates of GW from Populations of Compact Binaries - Spatial Distributions 



reference 



model 



^disk -Rdisk 23 char (disk) r bulge 
(kpc) (kpc) (kpc) (kpc) 



this work, Nelemans et al. ( 2004 1 



see Sec. 4.1 



this work (thick disk) 



0.3 


2.5 1 


6.4 2 


0.5 


1.25 


2.5 1 


6.8 2 


0.5 


0.2 


2.5 


5.3 





Nelemans et al. (2001b) 
t 



Edlund et al.M2005T 



e -«/«di* sech 2 (Z/Z disk ) 



Timpano et al. ( 2006 ) 



see|Hils fc Bender ( 2000 i and Nelemans et al.|( 2004 



Seto 1 2002 I 



e —R/ H diske — z/z disk 



0.09 



3.5 



4.8 



Bender & Hils (|1997 1 
H 



Hils et al. dl990) 



e -R/ Rdisk c -z/z disk 009 35 4< , 
e -R/R diske -z/Z diBk o.09 3.5 4.1 



Hils & Bender ( 2000 ) 



Hils ( 1998 i 



us 
us 



11S 
11S 



IIS 

ns 



us 
us 



Yu & Jeffery 1 2010 I 



e -fl/fldi,k sech 2 (Z/Z disk ) 0.352 2.5 5.6 
c -R/R diskc -z/z disk 1A58 2 ^ 5 6 ^ 5 



0.5 



Liu et al. 1 2010 1; Liu 1 2009 I 



-R/ -Rdisk e -z/z di3k 



0.09 



2.5 



5.0 



Farmer & Phinney 1 2003 I 



Webbink & Han 1 1998 1 



e — R/ Rdisk e — Z/Z diBk 



0.09 4.0 



4.7 



Ruiter et al. 



Ruiter et al. 



(2010 



(2009 



e — R/ Rdisk e — Z/^disk 

(halo) 



0.2 



2.5 



5.4 



0.5 



Benacquista et al. ( 2004 i 



c --R/-Rd 1B k sech 2 (Z/Z diak ) 



0.2 


2.5 


5.3 




4.2 


5 


7.2 




0.46 


1.6 


5.7 


see paper 


1.2 


10 


9.6 





Postnov & Prokhorov ( 1998 ) 



Lipunov et al. (1995 



Lipunov & Postnov (1987) 



Lipunov et al. (1987) 



q — R/ Rdisk C — (Z/Z disk ) 
e ~R/ Hdisk e — ^/-^disk 

e -B/B di s ke -(Z/Zdisk) 2 



Giampieri & Polnarev dj 1997 1 
Evans et al. (1987) 



q — R/ Rdisk C — (Z/Z disk ) 



0.1 



3.5 



4.5 



constant density 



0.2 



10 



Notes: ns = not specified. Distance to the Galactic Center is assumed to be 8.5 kpc for all models. 
Exponential squared bulge unless otherwise specified: e - ' r / rbul s°) with X> c har (bulge) = 8.4kpc. 
1 Numerical fit, DDWD systems only. 

bulge for this case. 



1. Assumptions about the relative numbers of systems in terms of total mass, mass ratio, and separation. These 
arise from different treatments of compact-binary progenitor evolution and of the stability of mass-transferring 
systems. Approaches include both analytic estimates and Monte Carlo population-synthesis calculations. 

2. Assumptions about Galactic structure and population space densities. These can affect estimates in two ways: 
first, the magnitude of GW signals will depend on the scale size assumed for the Galactic population; second, the 
"calibration" of the local number density of sources using observations of a particular source type will depend 
sensitively on the assumptions about Galactic structure. 

We now present a semi-analytical derivation of the amplitude and s hape of the GW continuum. Similar model s have 
been used extensively in GW data analysis, see for example works by Hughes ( 2002 ) ; Barack & Cutler (2004a b). Our 
results show that the magnitude and spectral shape of the continuum are a robust function ot the merger rate and 
characteristic chirp mass of Galactic bi nary systems, and d epen d only weakly on Galactic str ucture. 

The GW lum inosity of 
b are given by (Peters & 



We follow the arguments presented in Evans et al 
a circular binar y and the change in its orbital energy as a 
Mathews|p63] ) 

dE 32 



(1987} and |Giampieri fc Polnare^ | (|1997[) 
rmction of GW frequency / 



dt 



10/3 jlO/3 



dE 

w 



5/3 f-1/3. 



(Al) 



where E is the energy of the system and / is the GW frequency. Combining these expressions yields the frequency 
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evolution of the system, 



'1 
dt 



^8/3^^5/3^1/3 



(A2) 



If 7£birth(/>r) is the birth rate per unit frequency and volume of systems with / or b < f/2, in a stationary regime the 
density of sources must be 



d 2 N 
djdx 



gjirthC/yr) _ 5 ^-8/3 c- M _ 5/3 



df/dt 



96 



a 5 



^birth(/,r)r n / 3 . 



(A3) 



By continuity, 7\Lbirth(/, r) must also be equal to the merger rate, 7?. mrg (/, r), of systems which merge at / or b > f/2. 
DDWD systems, which generate most of the GW foreground, merge at / > 10~ 2 Hz, and most of them are born below 
a few 10 -4 Hz (although estimates vary). Thus, if we restrict Eq. (A3) between those two frequencies, we may just 
set 7?-birth(/, r) = ^mrg(r), defined as the total merger rate above 10~ z Hz. 
Now, the direction-averaged GW strain from a circular binary is given by 



h - { h + + h x )- n2d2c3 I df 



32 



.4/3 



QW/3 



M 



10/3 



/4 /3 



(A4) 



c° 

see for instance (Giampieri & Polnarev 1997[ ), where d is the luminosity distance to the binary; the corresponding 
detector response is 

h 2 j = (h 2 + Fl + h 2 x F*) = h 2 F 2 (r,t), (A5) 



and 



where F + and F x are the antenna-pattern functions (see, e.g., Maggiore 2007) and F 2 (r,t) is the polarization 
inclination-average d antenn a pat tern ((FT) + (F 2 ))/2 in the direction f. 

We combine Eqs. A3 A4 and |A5| to obtain the one-sided detector response density Shj- The yearly orbital motion 
of LISA-like detectors modulates the amplitude of signals to individual so urces; this effect averages out for isotropic 
sourc e distributions, but not for the anisotropic Galactic-binary foreground ( |Giampieri fc Polnarev|1997 Edlund et al. 
2005 1 . We therefore average over a year (while assuming for now that all binaries have the same chirp mass M z ), to 
obtain 



S hl (/) = (/ h 2 (d(v)J) F 2 (r,t) 



— dr 
dfdv i t 



J 7r -4/3^ Mx 5/3 / -7/3 



(F 2 (v,t)) t TZ mrg (v) 
c J J d(r) 2 

l 7r -4/3^! Mx 5/3 / -7/3 ( ^ )f t( ^ mrg/d 2 )r _ 



dr 



(A6) 



Of the two quantities defined in the last line of Eq. (A6), (-F 2 )?,* accounts for the time-, merger-rate- and direction- 
averaged instrument response, and it depends on the orbit and geometry of the detector. For an isotropic distribution 
of sources and a triangular interferometer with 60° an gles, F 2 (r,t) = 1/5 x sin 2 (60 °) = 3/20, which is often taken as 
the average isotropic response for LISA-like detectors ( Giampieri fc Polnarev||1997[). The factor 1/5 is the average of 
((FVj + (F 2 ))/2 over source locations and polarizations (again, e.g., in |Maggiore|2007| . Again for a LISA-like detector, 
but for typical disk populations, F 2 (r,t) is found empirically to vary by ~ 10% from 3/20. 

The quantity (7?. mr g/d 2 ) r is the merger-rate-averaged inverse squared distance. For a Galactic disk, it depends on 
the scale height and radial scale length, and of course on the position of the detector within the disk. A typical disk 
density distribution considered in many papers is the squared-hyperbolic-secant disk model, 



(l") 



7 ^e- R / R ^sech 2 (Z/Z disk ), 



(A7) 



where lZ ga ,i is the total merger rate in the Galaxy and V ga \ is the volume of the Galaxy. We find empirically that 
(R- mlg /d 2 } r varies slowly for 2.5 kpc < i?disk < 10 kpc and 0.2 kpc < Zdisk < 1 kpc. Table [5] summarizes the wide 
range of Galactic models used in the literature. For most disk models the characteristic distance to Galactic sources 
2? char = (V r ga i/7?. ga i) _1 / 2 (7^ mr g/c? 2 )~ 1 / 2 is approximately 5 to 6 kpc, and 2? c har + i?disk is close to 8.5 kpc, the assumed 
distance to the Galactic center. 



Inserting numerical values appropriate for our simulations (cases 1-3) in Eq. ( A6 1 yields 



S hl (f) ~ 1.9 x 10- 44 (//Hz) 



-7/3 



Hz" 



(F 2 



2? cha 



11 



gal 



M 



z, char 



5/3 



(A8) 



^6.4 kpc J VO-OWyr/ V°- 35M 0. 

where A^.char is the characteristic chirp mass (M.t ' 3 )~ 5 / 3 , The power spectral density Sh(f) of the GW strain is 
obtained by dropping the factor (F 2 )f jt (which is close to 3/20 for LISA-like detector). Equation (A8) is in good 
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agreement with the results of our simulations. Note however that the spectra of Fig. [2] are plotted for the TDI 
observable X (rather than for strain h or strain response hi), so they must be compared with 



S ^(f) ~ 1 . n a 2 'SWCO' 
1 + 0.6 x l 



(A9) 



where x = 2nfL/c, and L is the detector armlength. 

Our analysis indicates that Sh T {f) and Sh(f) are relatively insensitive to the structure of the Galactic disk, consistent 
with the results of our simulations, and instead depend strongly on the merger rate TZ ga \ and the characteristic chirp 
mass M z , c har- By constrast, we expect that the local population density of sources, as well as the number of resolvable 
source as a function of their GW strain, will depend roughly linearly on the disk's scale height. 

The transition between the confused and resolved detection regimes is characterized by dN/df, which is de term ined 
by the evolution of the chirp frequency and by the absolute merger-rate density of sources. Integrating Eq. (A3 1 over 
r we find 



dN 



5 x 1(T 3 (//Hz) -11 / 3 Has -1 



K 



gal 



0.015/yr/ V 0.35M( 



xlia 



-5/3 



(A10) 



consistent with our simulations. 

This simple analytical derivation of the GW continuum shows that the major differences between the various esti- 
mates listed in Table [4] arise primarily from two sources: the total merger rate of systems 1Z sa \ and the composition 
of the source population, leading to a different characteristic chirp mass Af ZiC h ar - £> c har is also relevant, but less 
important. (Remember however that our derivation is only predictive for frequencies high enough that most binaries 
form at longer periods — for LISA-like detectors, this is true above 10~ 4 Hz. For lower frequencies, the shape of the 
continuum depends on the detailed accounting of birth rate vs. period.) 



LOW-FREQUENCY GW INTERFEROMETERS AND THEIR RESPONSE TO BINARY WAVES 



LISA ([Prince et al.||2009| |Jennrich|[2009| , the space-based GW observatory studied jointly by NASA and ESA 
until early 2011, comprises three spacecraft orbiting the Sun in a quasi-equilateral configuration with armlengths 
~ 5 x 10 km. LISA measures GW strain by monitoring the oscillating distance between freely-falling test masses, 
using heterodyne laser interferometry. The spacecraft hover around the test masses to protect them from external 
perturbations, performing slight orbital corrections with extremely precise micro-Newton thrusters. 

Although LISA is sometimes described as a large Michelson interferometer in space, the lasers exchanged by space- 
craft are neither split nor reflected. Rather, two lasers on each spacecraft are used to establish six one-way links; each 
link provides a test-mass-referenced interferometric measurement y^ of the phase (or equivalently, the frequency) of 
the incoming laser, compared to the local, outgoing beam. In terms of fractional la ser frequency, the effect of a plane 
GW on the links is given by the two-pulse response ( Estabrook fc Wahlquist||l975 ) 



1 hjj ■ [h(pj,t~ Ljj) - h(pj,t)] ■ hi 

2 1 — hij ■ k 



(Bl) 



where nh is the unit vector from spacecraft i to j; k the GW propagation vector; h(pk,t) the TT-gauge GW strain 



(see, e.g., Maggiore 2007) at the spacecraft position p^ and at time t; and the photon time of flight along the 
arm. Thus, a GW pulse is registered twice in y^: once when it impinges on spacecraft j, and a time L+j after it has 
impinged on spacecraft i. In the limit of long GW wavelengths A ^> L, the response to plane GWs is therefore 
proportional to L. 

Each measurement is also affected by displacement noise from residual test-mass accelerations, position noise 
from optical-path and shot noise in the heterodyne measurements, and, most importantly, by oscillations in the 
central laser frequency, which are many orders of magnitude stronger than expected GW imprints, and which appear 
as y^ = Cjiit) — Cij(t — Lij), with CV,- the frequency noises for the six lasers. Laser frequency noise is removed by 
combining multiple one-way measurements with appropriate time delays so that each laser appears in canceling pairs, 
as in the unequal-arm Michelson observable 

L21)] 
Lsi)]. 



X(t) =[y 3 i(t) + y l3 (t 

-[«2l(t) + !/12(t 



L31) 
L21) 



2/21 (i 
2/31 (t 



L31 
L21 



Lis) 
L12) 



2/13 (* 



L31 
L21 



L13 

L\2 



(B2) 



This nois e-suppression te chnique, specific to LISA-like interferometers, is known as time-delay interferometry (TDI; 



see, e.g., 



Vallisneri 



2005 



for a review). TDI observables can be seen as synthesized interferometers, since their y^ 
components trace the paths (and reproduce the phase delays) of light across recognizable interferometric configurations 
(hence X's "Michelson" designation). Many different observables are possible with LISA's six laser links, but all can 
be reconstructed from a set of three o bservables, such as the unequal-arm Michelson X, Y, and Z centered on the 
three spacecraft ( Vallisneri et al.p008 ). 
In this paper, we assume an rms displacement noise 



Slg(f) = 3 x 10- 15 ms- 2 Hz 



-1/2 



/ 



10~ 4 Hz 



0.008 Hz 



(B3) 
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for each of the proof masses, and a total rms position noise 



«f(/, = 18*10-"„H z -^ 1 + (_L_) 



(B4) 



for each measurement. 

The GW mission currently under study as an ESA-led project ( Amaro-Seoane et al. 2012) would decrease cost 
primarily by reducing spacecraft and propellant mass, allowing for launch on smaller rockets than envisaged for LISA. 
Propellant is saved by placing the spacecraft closer together and closer to Earth. At the low-frequency end, the shorter 
armlengths reduce the response to GWs proportionally; for the same displacement noise, sensitivity also decreases by 
the same ratio. At the high-frequency end, the laser power available for position measurement increases as L~ 2 , since 
beams are broadly defocused at millions of kms, improving shot noise (but not other optical noises) by an rms factor 
As a consequence, the sweet spot of the LISA sensitivity shifts to higher frequencies. 

Spacecraft mass may also be saved by adopting a two-arm "mother-daughters" configuration with only four laser 
links, two each between the middle "mother" spacecraft and one of the "daughters". With four links, only one 
independent TDI observable can be formed, with a loss in sensitivity averaging \/2 at compact-binary frequencies. 
Because no trade studies have yet been performed to optimize the payload for the new orbits, in this paper we consider 
a "short LISA" configuration with aggressively reduced 1 x 10 6 -km armlengths, and with the LISA displacement and 
position noises (except for the power rescaling of the position-noise component, which we assume as five times smaller 
in rms). The actual sensitivity of a cost-constrained space-based GW observatory is likely to fall somewhere between 
this short LISA and the classic configuration. 

Searches for GWs of known shape proceed by correlating detector data with theoretical waveform templates computed 
for the expected range of source parameters. The statistical confidence of detection for a GW signal h immersed in 
Gaussian additive noise n is described by the signal-to-noise ratio (SNR) of the da ta s = h + n after filterin g by the 
best-fitting normalized template t. For a single TDI observable X, this is given by ( |Cutler fc Flanagan||l994[ ) 



p(s; i) = 4 Re 



Sx{f) 



df, 



(B5) 



where X s (f) and X$(f) are the Fourier transforms of the TDI data and of the TDI response to waveform t respectively; 
"*" denotes the complex conjugate; and Sx(f) is the noise power spectral density of TDI observable X. Normalization 
amounts to requiring p{t;t) = 1. 

With this noise weighting, the false-alarm probability of SNR exceeding a chosen detection threshold pt\ u - is pro- 
portional to exp(— pj hr /2). Because the many templates used in typical searches amount to repeated independent 
trials, pthr is routinely set to relatively large values (e.g., 5-10) to yield one false alarm per many years of data. The 
detectability of a signal h is characterized by its optimal SNR p(h;h), which is the SNR that would be obtained on 
average over noise realizations for a perfectly matching template. 
The noise Sx(f) includes the confusion- noise foreground Sjf nf (/) from unresolved sources (see Sec. [5]), as well 

For the unequal-arm Michelson observable X in the (rather reasonable) limit of equal 



as instrument noise Sx st 



interferometer armlengths, S^ st is given by ( Vallisneri et al.|[2008 ) 



S 



inst 
X 



In these units S, 



MIS 



16sin 2 (27rL/) 
MIS, 



2 1 



(2ttL/))S^(/) + 



(B6) 



= S pm (f)/(2wcf), So P J,J (/) = 5 op (/) x (2tt/ 



To obtain the TDI response to a GW signal Eq. (15j) , we i nsert Eq. (15j) , multiplied by the appropriate phasing function 
and sky-position-dependent polarization tensors, into Eq. ( |B2[ ). The response is modulated further by the heliocentric 
motion of the detector (which yields a time-dependent Doppler shift) and its yearly rotation in a plane inclined 60° 
from the ecliptic (which generates sidebands at frequency multiples of l/yr)Q Both effects imprint sky-position and 
polarization angles into the TDI response; this structure makes it possible to resolve sufficiently strong GW s ources 
that overlap in frequency space, and to determine their approximate sky position (see, e.g., Krolak et al.|[20 04). The 
TDI response to a population of sources is given simply by the linear superposition of the responses to individual 
sources. 



FITS OF THE PARTIALLY SUBTRACTED CONFUSION-NOISE FOREGROUND 



Barack and Cutler (2004a 2004b) codified the lore about the total noise of a LISA-like detector (i.e., the instrument 
noise plus Galactic-binary confusion noise) by providing a heuristic derivation for the well-known expression 



S] 



total.BC 



where they estimate 



(/) = min^a) x ^^W, S^(f) + S] 



£Conf,BC (/) = lA x 10-44(//H Z )~ 7/3 H Z -\ 



conf,BC 



(/)], 



(CI) 



(C2) 



7 Because "short LISA" orbits are very similar to those planned 
for classic LISA, this discussion is valid for both. Different orbits, 



such as halo trajectories around Lagrange points or Earth-centric 
orbits, would yield different signal modulations. 
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^ = 2xlO- 3 (//Hz)- 11/3 Hz-\ (C3) 

and kT^ 1 ~ 1.5/yr. As we discuss in App.JXJ under generic assumptions the unsubtracted Galactic-binary foreground 
will follow a /~ 7 / 3 spectrum; this is confirmed in our simulations, with a coefficient 2.6 x 1CP 44 for cases 1-3 and 5, 
and five times smaller for case 4 (see Fig.l3|. After the subtraction of detectable sources, the total noise found in our 
simulations is fit better by a smaller amplitude, a slightly different exponent, and a smoother transition function: 

£totaLsim (/) = ^inst + ^conLsim ^ tanh « ( 1^ dJV/d/) , (C4) 

with 

^conf,sim (/) = 14 x 10 -45( //Hz )-«/ 3 Hz -l ! ( C5 ) 



dA 



= 5 x 10" 3 (//Hz) 



~ 11/3 R Z -K 



(C6) 



For cases 1-3 and 5, and for our reference subtraction run (one year, one interferometric observable, pthr = 5), taking 
a = 1, j3 = 1/yr yields a good fit for a 5 Mkm LISA-like detector, while a = 0.7, f3 = 1.2/yr is appropriate for 1 Mkm 
arms. For the entire range of p c g that we probed, good fits are obtained by setting 



O 



2-1.2(p off /5)+0.2(p off /5) 2 , p 
1.6-1.2(p cff /5) + 0.3(p cfr /5) 2 , 



[-0.2 + 1.5 (p off /5) - 0.3 (p cff /5) 2 ]/yr 
[0.3 + 1.3 (p eff /5) - 0.4 (p off /5) 2 ]/yr 



(5 Mkm), 
(lMkm). 



(C7) 



For case 4, an acceptable fit follows from reducing gconf.sim by a factor of five. 



The equivalent confusion noise in the TDI X combination (expressed as fractional-frequency fluctuations) is given 



by 



ototal,sim/ r\ _ 3 16 £ sin X „total,; 

Uj "20 1 + 0.6 a; 2 h 



'(/), with x = 2ir fL/c. 



(C8) 



Here 3/20 is the all-sky averaging (F 2 )f yt discussed above; 1 + 0.6 x 2 is an approximation to the frequency-dependent 
interfer ometer response function; and 16 a: 2 sin 2 a; is the TDI X transfer function. [In Barack and Cutler's convention 
(2004b), which we follow here, the average response function is included in the definition of the equivalent strain noise 



5totai ) SQ that the s ky-averaged SNR can be computed directly as 4 Re J \h\ 2 / 'S£ otaJ (/) d/; therefore S£ onf , rather than 
Sf^ nl , appears in 5^ otal . By contrast, TDI observables already include the position-dependent instrument response, so 
it is appropriate to define the TDI noise S^ otal without the sky averaging.] 

Due to the evolving orientation of the LISA-like formation with respect to the Galaxy, the Galactic foreg round (before 
and after subtractio n) has a strong seasonal v ariation in amplitude, almost constant across frequencies ( Giampieri fc| 
Polnarev||1997| [Seto|[2004{ |Edlund et al.|2005[ ). For a single TDI X combination, for both our detector configurations, 
for all scenarios, and tor our assumpti ons "about the relativ e orientation of detector and Galaxy at t = [consistent 
with the standard MLDC LISA orbits ( Arnaud et al.| 2006)], this variation is fit well by 



-rconr.sim 



conr ,sim 



(/) x [l - 0.21cos(27rt/yr) - 0.45cos(47rt/yr) 



-0.05 sin(27r t/yr) + 0.01 sin(47r t/yr)] , 



(C9) 



where St° nt ' Slm (f , t) denotes the power spectral density computed for data segments s hort er than one year, but longer 



than the typical measurement timescale of a few hours. An equation analogous to Eq. (C9), with the same coefficients, 
holds also for the unsubtracted background Sh(f,t). 
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